Measurement of the atom number distribution in an optical tweezer using 

single photon counting 
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We demonstrate in this paper a method to reconstruct the atom number distribution of a cloud 
containing a few tens of cold atoms. The atoms are first loaded from a magneto-optical trap into 
a microscopic optical dipole trap and then released in a resonant light probe where they undergo 
a Brownian motion and scatter photons. We count the number of photon events detected on an 
image intensifier. Using the response of our detection system to a single atom as a calibration, we 
extract the atom number distribution when the trap is loaded with more than one atom. The atom 
number distribution is found to be compatible with a Poisson distribution. 
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I. INTRODUCTION 

The last few years have seen a growing interest in the 
study of mesoscopic systems consisting of typically a few 
tens of interacting particles. The properties of these sys- 
tems usually cannot be described by a mean-field ap- 
proach and are already too complicated to be calculated 
from the behavior of each individual interacting particle. 
Dense clouds of ultra-cold atoms provide an ideal test 
bed to study these mesoscopic systems. The interactions 
are well-understood at the two-body level and, experi- 
mentally, one benefits from the host of tools developed 
over the years to investigate ultra-cold atomic clouds (see 
IH)- Various implementations have already been re- 
alized such as arrays of optical tweezers 0] or of magnetic 
traps [s^l , double- well potential geometries [5] and optical 
lattices (e.g.[5i, i6j]). 

Beyond their intrinsic interest, a possible application 
of dense clouds of ultra-cold atoms is the production of 
atomic states for which the fluctuations of the number of 
atoms are reduced with respect to the Poissonian case. 
These clouds could for example be useful for atomic inter- 
ferometry below the standard quantum limit Recent 
demonstrations have been achieved using Bose-Einstein 
condensates with a few thousand atoms ([8j: SI)- The 
clouds could also be used as a source delivering a given 
number of atoms, as was already demonstrated in the 
single atom case in an optical tweezer ^lOj or in a one- 
dimensional optical lattice [ll|. For larger numbers of 
atoms, experiments showing the reduction of the atom- 
number fluctuations have bee n rep orted in Bose-Einstein 
condensates in an optical trap [l2| . in optical lattices [l^ , 
or in arrays of magnetic micro-traps [ij] . For all these ap- 
plications the knowledge of the distribution of the num- 
ber of atoms, and not only the average number of atoms, 
is essential. 

The number of atoms in a cold atomic cloud can be 
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measured by usual optical techniques, such as laser in- 
duced fluorescence or absorption. Nevertheless, special 
care must be taken when the number of interacting par- 
ticles involved is small. For instance, when using fluo- 
rescence imaging the small number of atoms combined 
with the low collection efficiency makes it hard to col- 
lect more than a few photons per realization of the ex- 
periment, therefore preventing the reliable extraction of 
the number of atoms in single shot. To circumvent this 
problem, a method was recently demonstrated [15] where 
freely propagating atoms fall through a sheet of resonant 
light leading to the detection of many photons per atom. 

In this paper we present a method to measure the dis- 
tribution of the number of atoms trapped in a micro- 
scopic optical dipole trap in a regime where the number 
of atoms is a few tens. The method relies on single- 
photon counting. The principle is as follows: we release 
the cloud of cold atoms from the dipole trap and let them 
expand in a resonant light probe. The atoms diffuse in 
the probe in a Brownian motion and emit photons. Some 
of the photons are collected by an imaging system and 
impinge on an image intensifier where we count them. 
By repeating the experiment several times and recording 
the counting results, we reconstruct the distribution of 
the number of detected photons. Our experimental setup 
can work either with exactly one atom, or with a cloud 
of up to a few tens of atoms. In this way we perform 
a calibration of the number of detected photon events 
when a single atom is trapped, and use this calibration 
to determine the atom number distribution. 

The paper is organized as follows. In Sect. |ll] we de- 
scribe briefly the setup. Section Hill describes the loading 
of our microscopic dipole trap in the single-atom regime 
and in the many-atom regime. In section ITVl we describe 
the experimental sequence and in Sect. |V] we present 
the image analysis used to count the number of photon 
events. Sect. IVII deals with the calibration of our detec- 
tion system using a single atom. Section|VlT]explains how 
we extract the atom number distribution in the multi- 
atom regime. In Sect. IVIIII we compare the mean num- 
ber of atoms extracted by this counting method to the 
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result obtained by direct integration of the laser induced 
fluorescence, and find good agreement between the two 
methods. 



II. EXPERIMENTAL SETUP 

Our experimental setup has been described in details 
elsewhere EMI and consists essentially in a microscopic 
optical dipole trap, in which we trap up to a few tens 
of ^^Rb atoms (see details in Section IIIII) . The trap is 
produced by a laser beam at 850 nm that is focused by 
an aspheric lens with a numerical aperture N.A. = 0.5. 
The resulting spot has a Gaussian profile with a waist 
w — l.f /im. A power of 9mW yields a trap depth 
of 2 mK. The dipole trap is loaded from a cold atomic 
cloud held in a magneto-optical trap produced around 
the focal point of the aspheric lens. The magneto-optical 
trap captures atoms from a Zeeman decelerated atomic 
beam coming from an oven. 

We use the same aspheric lens to collect the fluores- 
cence light at 780 nm. A dichroic beam-splitter sepa- 
rates the fluorescence light from the dipole trap light. A 
fraction of the fluorescence light is sent to a fiber-coupled 
avalanche photodiode operating in a single photon count- 
ing mode. The remaining fluorescence is directed onto an 
image intensifier, followed by a low-noise CCD camera. 
We have measured a detection efficiency rj ~ 3 x 10~^, 
which takes into account the solid angle of the aspheric 
lens (7%), the transmission of the optics (40%), and the 
measured quantum efficiency of the intensifier photocath- 
ode (10%). Here, we note that using an image intensi- 
fier is crucial for the experiment presented in this paper 
because it enables our imaging system to detect single 
photon events. In the absence of the intensifier, single 
photon events have an amplitude lower than the noise of 
the CCD camera alone (6e~ /pixel). The intensifier am- 
plifies a single photon event to an average amplitude of 
920e~ /pixel, well above this noise level. 

In order to characterize the atomic cloud, we use a 
probe consisting in two counter-propagating laser beams 
in a (7_|_ — (T_ configuration. This probe is tuned to 
the {5'^Si/2,F = 2) to {5^P3/2,F' = 3) transition and 
is sent together with a repumping light tuned to the 
{5'^Si/2,F = 1) to (52P3/2, F' = 2) transition. The probe 
light is focused down to a waist of 0.9 mm and has a sat- 
uration parameter ///sat ~ 2.7 for each beam. 



the loading from the magneto-optical trap, and the two- 
body losses due to collisions assisted by the near-resonant 
light of the cooling and trapping lasers [l^. When the 
loading rate of the dipole trap is lower than the light- 
assisted collision rate, only one atom can be present at a 
time in the trap: if a second atom enters the trap a light- 
assisted collision takes place and both atoms are lost. On 
the contrary, when the loading rate of the dipole trap 
overcomes the inelastic loss rate, the average number of 
atoms is larger than one [l^, [lo] ■ For the experiments 
presented below we find that, after a fixed loading time 
(;$ 1 s typically), the mean number of atoms in our mi- 
croscopic trap increases linearly with the trap depth and 
reaches N — 9 for a trap depth of ^ 20 mK. We mea- 
sure the temperature of the trapped atom(s) using either 
a release- and-recapture method [17] or a time-of-fiight 
method [l^. We find a temperature around 150 /iK in 
the single atom case and around 1 mK in the multi-atom 
regime. 

In the single-atom regime, the presence of an atom in 
the trap is revealed b y a step in the fluorescence signal 
collected on the APD [Ig . This step is used to trigger an 
imaging sequence where the number of fluorescing atoms 
is exactly 1: this way there are no fluctuations in the 
atom number from shot to shot (variance AA^^ = 0). In 
the multi-atom regime, the situation is different as we 
observe no such steps on the APD signal: firstly because 
when the atomic density increases, the time spent on av- 
erage by each atom in the trap decreases below the time 
bin resolution of the APD counter due to an increasing 
inelastic loss rate, and secondly because the height of the 
steps itself decreases due to light-induced dipole-dipole 
interactions that shift and broaden the atomic line 2^]. 
One thus cannot rely on fiuorescence steps to perform an 
in situ measurement of the number of trapped atoms. In- 
stead, the technique presented in this paper goes around 
this problem by converting the multi-atom signal into a 
sum of single atom events, spatially separated on a CCD 
camera. In this second regime, the absence of fluores- 
cence steps on the APD signal also makes it impossible 
to trigger the imaging sequence after a given number of 
atoms is loaded in the trap. We thus trigger the imag- 
ing sequence after a flxed loading duration. Due to the 
random nature of the two-body losses during the load- 
ing this procedure leads to fluctuations in the number of 
trapped atoms from shot to shot (AA^^ ^ 0) [l9| . 



III. LOADING OF THE DIPOLE TRAP 



IV. EXPERIMENTAL PROCEDURE AND 
FLUORESCENCE IMAGES 



The experiment can operate either in a regime where 
only one atom is present in the dipole trap or in a regime 
where up to 10 atoms are trapped. The transition be- 
tween the two regimes is governed by the atomic density 
in the magneto-optical trap, which is controlled by the 
flux of the atomic beam. The average number of atoms 
N in the dipole trap is set by two competing processes: 



Once the trap is loaded, either with exactly one atom, 
or with N atoms, we switch off the cooling and repump- 
ing lasers and wait an additional 30 ms before switching 
off the dipole trap. This waiting period allows enough 
time for the atoms of the magneto-optical trap to leave 
the observation region, otherwise they would strongly 
contribute to the signal observed with the intensifier. We 
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FIG. 1: (Color Online) Figure (a) is an image taken for a 40^s- 
illumination period of a single atom. The region of interest 
has a size of 50 /^mx50 /im in the plane of the atoms, corre- 
sponding to a measured temperature of 150 fiK (see text). 
In figure (b), about 8 atoms were released from the trap 
and illuminated with a 20 /is probe pulse. The image size 
is 100 /^mxlOO ^m in the plane of the atoms, corresponding 
to a temperature of 1 mK. 



then switch off the dipole trap (in ~ 200 ns) and, after a 
time of flight of the atoms of 1 fis in free space, we turn on 
the resonant probe, while gating the intensifier for a given 
period of time At (in the 10 — 60 /isec range). Finally, we 
read out the CCD camera, which lasts ~ 500 ms, while 
launching the next loading sequence. 

During the 1 /is time of flight in free space, the atoms 
start separating from each other. Then they walk ran- 
domly in the retro-reflected light probe and scatter pho- 
tons at a rate ^ T/2 ~ 2 x 10'' s^^ (F is the line-width 
of the atomic transition) . Some of the scattered photons 
are collected on the image intensifier and impinge on the 
CCD at different positions due to the random walk of 
the atom in the probe light [22i]. The total number n of 
detected events is proportional to the number of atoms 
N on average. 

The images shown in Fig. [ija) and (b) were obtained 
after a single realization of the experimental sequence. 
Figure [TJa) corresponds to the time of flight of a single 
atom moving randomly in the probe light. Figure [Hb) 
corresponds to the time of flight of about 8 initially 
trapped atoms. In both cases, the detected photons 
correspond in a large extent to photons being scattered 
by the atom and detected by the intensifler and, in a 
smaller extent, to photons due to spurious laser light 
or self-induced charges generated inside the intensifier 
(see more details below). This fact led us to adjust the 
size of the region of interest to the measured temper- 
ature of the atom(s): for a time of flight of duration 
At, we chose a square with a side dimension ^ 4(T^Ai, 
where — yJk-QT/M, T being the temperature of the 
atoms, M their mass, and the Boltzman constant. 
This choice leads to a negligible probability of missing an 
atomic event, while reducing the number of background 
events, which are uniformly distributed on the CCD. 



V. PROCEDURE FOR ANALYZING THE 
IMAGES 

In order to count the number of events detected on the 
CCD camera, one has to decide what one calls an event. 
In order to do so, we flrst characterize the noise of the 
CCD camera, which is essentially due to the readout of 
the CCD chip. We flnd a root-mean-square (rms) noise 
fccD — 6e^ /pixel. We then apply the following proce- 
dure. We search for the pixel with the largest amplitude 
of the signal. If the signal amplitude on this pixel is 
larger than 6crccD, we call it a peak and exclude from 
the image a zone of 13 x 13 pixels around it, the rms size 
of an individual event being 1.3 x 1.3 pixels (set by the 
intensifler resolution and corresponding to 0.6 x 0.6 /xm^ 
in the plane of the atoms) [1^ . We then repeat this pro- 
cedure by searching for the pixel with the next largest 
amplitude until no peak with an amplitude larger than 
6crccD is found. 

We have independently measured the amplitude dis- 
tribution of single-photon events recorded by the inten- 
sifler followed by the CCD camera. Using this distribu- 
tion, we have found that less than 4% of the events have 
a peak amplitude smaller than 6(Tccd- Our procedure 
for counting the events therefore does not underestimate 
their number by more than 4%. 



VI. CALIBRATION USING A SINGLE ATOM 

The flrst step of the calibration of the imaging sys- 
tem as a counting device was to measure the histogram 
of the number of background events in one image when 
no atom is trapped in the dipole trap. The background 
events have three different origins: the largest contri- 
bution comes from light scattering by the atoms of the 
atomic beam that cross the region of observation (60%); 
the second largest contribution comes from self-induced 
events generated by the gated intensifler when the high 
voltage is switched on and off (33% of the counts); flnally, 
a smaller contribution comes from scattering of the probe 
beams on the surfaces inside the vacuum chamber (7%). 
The histogram of the total number of background events 
is shown in Fig. [51 which results from the analysis of 
200 images of the background. The probe pulse duration 
was At — 30 /IS. The data are well fltted by a Poisson 
distribution with average nbg = 0.7. 

In a second step, we repeated the same experiment 
with exactly one atom trapped in the dipole trap. In 
particular, the probe duration was the same as above 
{At = 30 /is). As the detection efficiency is small and 
corresponds to a random collection of photons, we also 
expect a Poisson distribution for the number of detected 
photons emitted by the atom. In this case, the distribu- 
tion of detected events, including the background events, 
is a Poisson distribution with a mean value: 

nsingle = ^at + ^bg ■ (1) 
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FIG. 2: (Color Online) Probability distribution of the num- 
ber of detected events in one image. The probe duration is 
30 Red diamonds : background events (no atom present 
in the trap). Blue circles : total number of detected events 
(including background events) when the dipole trap is filled 
with exactly one atom. In both cases, the probability p to 
measure a given number a counts in one image is deduced by 
analyzing 200 images. The error bar on each data point is 
calculated by — p)/200, which we tested by repeatabil- 

ity measurements. A Poisson distribution fits the data very 
well in both cases (dotted lines). 



The experimental distribution is shown in Fig. [51 to- 
gether with a fit by a Poisson distribution. We obtain 
a good agreement for fisingio = 2.5. This yields an aver- 
age number of detected photons, emitted by one atom, 
of nat = 1-8. 

Finally, we tested the linearity of our counting system 
with the duration of the probe pulse. The results are 
shown in Fig. |21 where we plot the average number of 
detected events due to scattering by a single atom, fiat, 
obtained after subtraction of the average number of back- 
ground events. We first performed the experiment with 
a single atom with a temperature T ~ 150 /iK and ob- 
served a linear dependency of fiat with the duration of 
the probe illumination. A linear fit to the data yields 
a photon detection rate, for one atom released from the 
trap, of i?d = 57200 ± 3100 s~^ (the error bar is from 
the fit). Knowing the detection efficiency of the imag- 
ing system ('--^ 3 x 10"'^), we deduce a scattering rate of 
1.9 X 10^ s^^, in agreement with the value calculated for 
the parameters of our probe. 

In order to further test the linearity of our counting 
system, we also performed the same measurement for 
an atom adiabatically cooled down to 16 /j,K [13] . In 
this case, the detected events tend to accumulate on a 
smaller area on the I-CCD, due to the smaller velocity 
of the atom (the time of flight of the atom is maintained 
constant and equal to 1 /iS with respect to the previous 
case at T ~ 150 /^K). In spite of the increase in the sur- 
face density of detected events, we were still able to dis- 
criminate between individual events and found the same 



FIG. 3: (Color Online) Average number of detected photons 
after scattering by a single atom with temperature 16 /iK 
(blue circles) or 150 /iK (red diamonds), versus the duration 
of the atom illumination by the probe light. Linear fits to 
the data (dashed lines) yield a photon detection rate _Rd = 
57200 ± 3100 s"^ for 16 /iK and i?d = 61800 ± 2400 s~^ for 
150 /iK. 



linear dependency as above, with a photon detection rate 
of 61800 ±2400s-i. 

To conclude, we detect on average fiat = 1.1 pho- 
tons scattered by a single atom for a probe duration 
Ai = 20 /iS. This value will serve as a calibration in the 
experiments described below, where the average number 
of atoms is the unknown. 



VII. ATOM NUMBER DISTRIBUTION IN THE 
MULTI-ATOM REGIME 

We now operate in the regime where more than one 
atom are loaded in the dipole trap. We repeat the 
same experimental procedure and extract the histogram 
of the number of detected events for a probe duration 
Ai = 20 /(S. Figure |4] shows the histogram of the back- 
ground events as well as the probability distribution of 
the total number of detected events (including the back- 
ground events) when about 6.3 atoms are loaded in the 
dipole trap on average. The background gives an aver- 
age number of events ffbg = 2.1, larger than in Sect. IVI[ 
as we needed to increase the flux of the atomic beam in 
order to load more than one atom in the trap. 

To fit the distribution of events in the multi-atom 
regime, we consider the distribution of the number of 
atoms N in the dipole trap, P{N,N), with N the mean 
atom number. The probability to detect n photon events 
is then given by the composed law p^ : 

oo 

p(n) = ^ P(iV, iV) X n(n, iV X fiat + fibg) (2) 

N=0 

with n(n, a) the Poisson distribution of mean a. Even in 
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FIG. 4: (Color Online) Probability distribution of the num- 
ber of detected events when the dipole trap is filled with 6.3 
atoms on average (blue circles). The data are well-fitted by 
a composed law that includes a Poisson law for the number 
of atoms in the trap. The data in red diamonds correspond 
to the background events only (fibg = 2.1) fitted by a Pois- 
son law. The duration of the probe is set to 20 fj,s so that 
the average number of detected photons scattered per atom 
is fiat = 1.1. 



the case when P{N, N) is a Poisson law, the distribution 
p{n) is not Poissonian. Nevertheless, for any distribution 
P{N,N), the mean value of this composed law is given 
by: 



^muiti = N X n^t+ ribg . 



(3) 



A direct calculation of the mean of the data shown in 
Fig. m yields fimuiti = 9. Taking into account that the 
probe duration has been chosen to detect fiat = 1.1 events 
per atom and that fibg = 2.1, equation ([3]) yields an av- 
erage number of atoms N = 6.3. Taking a Poisson dis- 
tribution for P{N, TV), a fit of the data by the composed 
law ([2]) leads to the same result. The result of the fit is 
shown in FigU) 

We now discuss the Poissonian assumption mentioned 
above for the atom number distribution P{N,N). We 
extract from the data the varian ce of the number of de- 
tected events, Anl^^^^,. = n^^^; - fi^^n^. This variance 
is related to the variance of the number of atoms AiV^ 
by the following expression, calculated using the proba- 
bility p{n) of equation ([2]) , and valid for any distribution 
P{N, N): 



X AN^ 



77-niulti 



(4) 



Taking again riat = 1.1, we find for the data of Fig. |4] 
a ratio of AN'^/N = 0.86 ± 0.13 (the error bar is sta- 
tistical). This value is compatible with P{N,N) being 
a Poisson distribution. We repeat this measurement for 
various average numbers of atoms ranging from TV = 2 to 
N = 9 and find that AN'^/N is equal to 0.76 in average 
with a rms dispersion of 0.13. This 1 a uncertainty does 
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FIG. 5: (Color Online) Comparison of the results for the av- 
erage number of atoms iV by two methods : direct counting 
of single fluorescence events (vertical axis) versus fluorescence 
integration (horizontal axis). Both methods give comparable 
results, as indicated by the linear flt to the data (dashed line) 
with a slope of ~ 1 (see text). The experiments were per- 
formed in two configurations: with a microscopic dipole trap 
(red diamonds, w = 1 fim) or with a larger dipole trap (blue 
dots, waist w — 4 fim). 



not clearly indicate a sub- Poissonian behavior as AN"^ /N 
does not deviate from 1 by more than 2 a. 



VIII. COMPARISON WITH FLUORESCENCE 
INTEGRATION 

We have checked that the average number of atoms ex- 
tracted by the counting method described so far is con- 
sistent with the number of atoms extracted by a second 
method based on fluorescence integration. This second 
method consists in illuminating the cloud of freely propa- 
gating atoms by a 2 /zs probe pulse, and accumulating the 
fluorescence detected on the CCD over many realizations 
of the experiment. In each time-of-flight experiment, we 
let the atoms fly during 1 /zs before we illuminate them 
with the probe pulse. As explained in [l^, after typically 
several hundreds repetitions of the experiment, the image 
shows a nearly gaussian distribution, reconstructed from 
many individual detected events. We perform this exper- 
iment in the multi-atom regime and then in the single- 
atom regime. We use the latter to extract the number of 
atoms in the multi-atom regime, by calculating the ratio 
of the integrated fluorescence values obtained in the two 
regimes. 

This fluorescence integration based method was imple- 
mented in two trap conflgurations that we now describe. 
As already explained, a significant fraction of the back- 
ground events come from photons being scattered when 
the atomic beam interacts with the probe light. To min- 
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imize this effect, we decreased the flux of the beam by 
reducing the temperature of the oven. This is at the ex- 
pense of the loading rate of the dipole trap: in this first 
configuration, we could not load more than '^10 atoms 
in the microscopic dipole trap. In order to calibrate our 
detection scheme with a larger number of atoms, we im- 
plemented a dipole trap with a larger waist w — 4: /xm. In 
this second configuration, the atomic density is smaller 
and the light-assisted collisions limit the number of atoms 
to a larger value than in the microscopic dipole trap [13] , 
typically 30. The results of the experiment are shown 
in Fig. [S] where we plot the average number of atoms 
obtained from the counting method versus the average 
number of atoms obtained by fluorescence integration. 
The numbers of atoms obtained by the two methods are 
compatible : a linear fit to the data yields a slope of 1.08 
which is compatible with 1, as the statistical uncertainty 
on the slope is 0.05 and we evaluate the bias uncertainties 
attached to each method to 4% for the counting method, 
and 11% for the integrated-fluorescence based method. 

IX. CONCLUSION 

We have demonstrated in this paper a method to 
reconstruct the distribution of the number of atoms 



trapped in an optical dipole trap. This method gives 
access to the distribution of the number of atoms for 
clouds containing up to a few tens of ultra-cold atoms. 
It is based on the detection, at the single photon level, of 
fiuorescence events scattered by the atoms when they un- 
dergo a random walk in the resonant probe light. So far, 
we have investigated a regime where the distribution of 
the number of atoms does not clearly depart from a Pois- 
son law. Our technique is in principle applicable to any 
distribution and could therefore be useful to investigate 
blockade effects leading to sub-Poissonian distributions 
of the number of trapped atoms. 
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